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Abstract. In a previous paper (arXiv: 1309.2204 JHEP 1311 (2013) 087), we present a holographic s+p 
superconductor model with a scalar triplet charged under an SU(2) gauge field in the bulk. We also study 
the competition and coexistence of the s-wave and p-wave orders in the probe limit. In this work we continue 
to study the model by considering the full back-reaction The model shows a rich phase structure and various 
condensate behaviors such as the “n-type” and “u-type” ones, which are also known as reentrant phase 
transitions in condensed matter physics. The phase transitions to the p-wave phase or s+p coexisting phase 
become first order in strong back-reaction cases. In these first order phase transitions, the free energy curve 
always forms a swallow tail shape, in which the unstable s+p solution can also play an important role. The 
phase diagrams of this model are given in terms of the dimension of the scalar order and the temperature 
in the cases of eight different values of the back reaction parameter, which show that the region for the 
s+p coexisting phase is enlarged with a small or medium back reaction parameter, but is reduced in the 
strong back-reaction cases. 


1 Introduction 

Over the past years,the AdS/CFT correspondence [HEE] 
has attracted a lot of attention as a tool to study strongly 
coupled systems. One of the most successful applications 
of the correspondence is the so-called holographic super¬ 
conductor which is first proposed in Refs. BU3- In this 
model, the condensation of a charged field in asymptoti¬ 
cally AdS spacetime can trigger the superconducting phase 
transition in the dual held theory. According to the dif¬ 
ferent property, especially the symmetry of the condensed 
fields, the holographic model has been extended to the 
p-wave and d-wave superconductor models Erasumn]. 

Recently the competition and coexistence between dif¬ 
ferent order parameters have been studied in some holo¬ 
graphic framework. In Refs. mmm, the authors stud¬ 
ied the system with two s-wave orders. Refs. mmm 
presented the results of the holographic systems with s- 
wave and p-wave orders simultaneously. The holographic 
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systems with one s-wave order and one d-wave order have 
also been investigated in Refs. [mug] , More works on 
holographic study of the multi-order systems can be found 
in Refs. These stud¬ 

ies might shed light to the understanding of some univer¬ 
sal properties of the competition and coexistence between 
different orders in real materials such as high T c cuprates 
and 3 He superfluid. 

In a previous work |14j . we built a holographic model 
with a scalar triplet charged under an SU(2) gauge held 
in the bulk and studied the competition and coexistence 
between the s-wave and p-wave superconducting order pa¬ 
rameters. The results showed that the system could be in 
an s-wave phase at higher temperatures and then trans¬ 
fers to a p-wave phase at lower temperatures via an s+p 
coexisting phase. In that work, we worked in the probe 
limit and found that the temperature region for the s+p 
coexisting phase is very narrow. The free energy curves of 
these phases show that the s+p phase is always stable. 

Our previous study in the probe limit has revealed a 
typical phase transition behavior among the s-wave phase, 
the p-wave phase and the s+p phase. However, the results 
in the probe limit can only be meaningful when the charge 
of the s-wave and p-wave orders is vary large. If one wants 
to uncover the complete phase diagram of this s+p model 
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with finite charge, we must go beyond the probe limit 
and take into account the back reaction of the matter 
fields on the metric. Moreover, note that in previous stud¬ 
ies, the back reacted s+s model has already shown new 
phase transition behaviors which do not exist in the probe 
limit PUDS- In the s+p model, the condensate of p-wave 
order would break the spatial rotational symmetry spon¬ 
taneously, which is obviously different from the s+s case. 
Thus we can expect to discover new phase transition be¬ 
haviors in the back reacted s+p model with an anisotropic 
p-wave order. Finally, the holographic p-wave supercon¬ 
ductors with back-reaction have already been well stud¬ 
ied it was found that the phase transition 

controlled by the p-wave order becomes first order when 
the back reaction is strong enough |311l32ll33] . Therefore, 
it is necessary to study whether there is still an s+p phase 
and what would happen to the s+p phase when the order 
of phase transition changes. In summary, it is interesting 
and important to take into account back-reaction in the 
s+p model to get the rich phenomena of phase transitions, 
and to reach a better understanding on systems with com¬ 
peting orders. 


In this paper, we turn on the back reaction of mat¬ 
ter fields on the background geometry in our holographic 
s+p model M- More precisely, the strength of the back 
reaction is controlled by a parameter b. For simplicity, we 
consider eight values of back reaction parameter b varying 
from 0.1 to 0.8, and build corresponding phase diagrams in 
terms of the dimension of the scalar order and the temper¬ 
ature. We also give an alternative setup of the holographic 
s+p system, which would allow us to vary the charge and 
dimension of the p-wave order. 


This article is organized as follows. In the next sec¬ 
tion, we briefly introduce our holographic model and the 
definition for the parameters. In section [3] we calculate 
the free energy of the system. In section we show the 
condensation behavior of some interesting phase transi¬ 
tions with different back reaction strength. We also give 
the free energy behavior for some special cases with unsta¬ 
ble s+p solution to illustrate the important role played by 
the unstable s+p solution. In section [5j we construct the 
phase diagrams with eight different values of back reaction 
parameter. In section [6j we give another setup of a holo¬ 
graphic s+p model, from which with some special choice of 
model parameters, we can get the same equations of mo¬ 
tions and thus the same results of phase transitions. We 
summarize our results and make some outlooks in section 

m 


2 Holographic model of an s+p 
superconductor 


The action of our holographic s+p superconductor model 
takes the following form [14] 

S = S G + S M , ( 1 ) 

Sg = J d d+1 x^g{R - 2T), (2) 

S M = 4 [ d d+1 xV=g(-D^ a D^ a 
9 C J 

~^F“„F apv - m 2 & a & a ). ( 3 ) 

Here F a is an SU(2) scalar triplet and is the SU( 2) 
gauge field with F“ the field strength, a = (1, 2,3) is the 
index of the generators of SU(2) algebra. From this action, 
we can get the equations of motion for the matter fields 

g^DpD„<F a - m 2 F a = 0, (4) 

g pv D p F^ - 2e abc 'F b D^ c = 0, (5) 

and for the gravitational fields 

R, iv - ^ (R~ 2A )9^ = ( 6 ) 

where b = n g /g c characterizes the strength of the back 
reaction of the matter field on the background geometry, 
and Tpv is the stress-energy tensor of the matter sector 

T,,u = (-D.rD" F a - \F“ v F a ^ - m 2 F a F a )g tlv 

-2 (~D^ a D^ a - \F; p F-P). (7) 

Here Dp is the covariant derivative with SU(2) connection 
on the curved spacetime background and e abc is antisym¬ 
metric with e 123 = 1. The action of Dp on a spacetime 
tensor T pi with one group index a is given by 


U T a 


— X7T a 
~ v m+i. 


+ £' 


ibc a b 


A°T C 


( 8 ) 


We will work in d = 3 case as in Ref. m- In order to 
get both the s-wave and p-wave superconducting phases 
in this model, we can choose the A p field as the electro¬ 
magnetic U{ 1) gauge field, and <F 3 , H 3 as the scalar and 
vector order parameters, respectively. Thus the ansatz for 
the matter fields is 


<F 3 = <F 3 (r), A 1 , = 0(r), A 3 = F x (r), (9) 


and all the other components of the scalar and gauge field 
are set to be zero. 

The black hole geometry compatible with the ansatz 
of matter fields @ is P1IMIIM] 

ds 2 = - N(r)a(r) 2 dt 2 + j^pjdr 2 

+r 2 f (r)~ 2 dx 2 + r 2 /(r) 2 dy 2 , 


( 10 ) 
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with 


r“ _ 2M(r) 


N(r) = ^(l- r 


), 


(ID 


where L is the AdS radius. The function M(r) is related 
to the mass and charge of the black brane. The horizon of 
the black brane is located at r = where M(rh ) = r b /2, 
and the temperature is 


T = 


N'(r h )a{r h ) 
47r 


( 12 ) 


The full equations of motion for this system under above 
ansatz read 

«'w = W! S^ {r? + j 

1 .2 r 2 2 2 r m 2 & 2 L 2 r 2 >F 3 (r) 2 </>(r) 2 

+-iww + 2N(r) l{ rf ' 

1 2 2 2 2 b 2 L 2 r 2 (b' (r) 2 

+ -bWiV(r)^(r) 2 + ^ 


+ 


cr'(r) = 


r 2 L 2 N(r)f'(r) 2 

2/(r) 2 ’ 

b 2 f(r) 2 ^ x {r) 2 (t)(r ) 2 5 2 /(r) 2 cr(r)iP r '(r) 

2rN(r) 2 a(r) 2 r 

6 2 r^ 3 (r) 2 </>(r) 2 2 , 2 

4- „ r/ NO / ' -h b ra{r)W 3 (r) 


(13) 


+ 


N(r) 2 a(r) 
rcr(r)/'(r) 2 

/W 2 ’ 


(14) 


= b 2 f{rf^ x {r) 2 (l){r) 2 + b 2 f (r) 3 V' x (r) 2 + f(r) 2 


2r 2 N(r) 2 a(r) 2 2 r 2 

f'(r)N'(r) f'(r)a'(r) 2/'(r) 


N(r) 


r(r) 


j.z /7 f(r) 2 ^ x {r) 2 (l)(r) 2<F 3 (r) 2 </(r) 

* (r > = —^- + JV(e 


+ 

*£(r) = - 


^"(r) = 


a'(r)(f>'(r) 2c)/(r) 

a(r) r ’ 

2f'(r)K(r) N'{r)iF x (r) 
f(r) N(r ) 

&x(r)(l>(r) 2 _ cr'(r)^(r) 
N(r) 2 a{r) 2 cr(r) ’ 

TO 2 >F 3 (r) ./V'(r)iF 3 (r) <F 3 (r)</(r) 2 

Af(r) A^(7’) N(r) 2 cr(r) 2 

cr'(r)^(r) _ 2^(r) 
cr(r) r 


/M 

(15) 


(16) 


(17) 


(18) 


Note that in above equations, there exist four sets of scal¬ 
ing symmetries as follows 

(1) TV A 2 7V , L -»■ A _1 L , b A _1 6 , a -»• A-V , 



(j) —i► X(f> , d'x —> \']/ x , 

'As —t AiF 3 , to —>■ Xm ; 

(19) 

(2) 

cf) — y A (f> , ^f x — y X&x i N — y N , 



M -> A 3 M , r — y Xr ; 

(20) 

(3) 

(j) — y A (j) , 

a —>• Act ; 

(21) 

(4) 



(22) 


These scaling symmetries are very useful in the numerical 
calculations. For example, we can use the first scaling sym¬ 
metry to set L = 1 and the second one to set = 1. After 
we get the numerical solutions, we can use these scaling 
symmetries again to recover L and to any value. The 
last two scaling symmetries will be used to scale the solu¬ 
tions to asymptotically AdS type with lim <r(r) —¥ 1 and 

r—too 

lim f(r) —> 1. 

r—±oo 

In order to solve the equations of motion numerically, 
we need to specify the boundary conditions both on the 
horizon r = and on the boundary r = oo. Without loss 
of generality, we set L = 1 in the rest of the paper. Then 
near the horizon, the functions can be expanded as 


M(r) = T f+M hl (r-r h ) + ..., (23) 

<r(r) = (T h o +<Thi(r ~ r h ) + ... , (24) 

f(r) = fho + fhi(r - r h ) + ... , (25) 

<t>(r) = (t’hiir - r h ) + (j) h2 (r - r h ) 2 + ... , (26) 

!?,(r) = ^xhO 4” ^xhlir t'h) 4“ , (27) 

&3(r) = &3h0 + ^3hi(r - r h ) + ... . (28) 


One can check that only the coefficients {a ho, fho, 4 > hi,'&xh0i ^ 3 / 10 } 
are independent. That means when we work at some fixed 
value of b and to 2 , we can get a group of solutions with 
these five independent parameters. However, we should 
further consider some constraints from the boundary r = 
oo. Then the degrees of freedom of the solutions would be 
reduced. 

We can expand the functions near the AdS boundary 
as 


M(r) 

= M b0 + ^ + .. 

f V ( 

■ , cr(r) = cr b0 + 

fir) 

= fbO 4- yf- + ... 

, t p 3 (r) = 

& x {r) 

= 'PxbO 4- 7 1 + ■■ 

• j = V- ~ 


0&3 

r 3 


&3E 

1 


(29) 


where A = (3 + y/9 + 4?n 2 )/2 is the scaling dimension of 
the scalar order. Note that there is no mass term for the 
SU(2) fields, thus the dimension for the p-wave order is 
fixed to be A p = d — 1 = 2. 

In order to make the boundary geometry to be asymp¬ 
totically AdS (a-AdS), we should have a b o = fbo = 1- 
These two conditions could be easily satisfied by a scale 
transformation from any known solution by using the last 
two scaling symmetries ( 21|22 ). 

In the asymptotically AdS spacetime, we can use the 
AdS/CFT correspondence to get information of the dual 
field theory. The AdS/CFT dictionary tells us that p and 
p are related to the chemical potential and charge density 
respectively, while <F 3 s, 'I , x b o are related to the sources and 
!F 3 e, &xbi are related to the expectation values of the dual 
s-wave and p-wave operators. To study the phases where 
the symmetry is spontaneously broken, we further impose 
the source free condition with i^o = ^ 3 s = 0. There is an 
alternative quantization with the expression of the source 
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and expectation interpretation exchanged 
But here we only focus on the standard quantization case. 

With these conditions from the boundary, we can count 
the degrees of freedom left for the equations. Besides the 
two parameters b and A from the model, we have five free 
parameters at horizon to fix the solution. The solution 
we need should satisfy four conditions on the boundary. 
Thus the solutions satisfying the four boundary conditions 
would have only one free parameter besides b and A. In 
this paper, we would choose this free parameter as T while 
fixing the value of /r. Thus we work in the grand canoni¬ 
cal ensemble, and we can plot our numerical results with 
respect to the dimensionless parameter T/p,. This is also 
convenient for comparing the results in this paper to that 
in the probe limit in Ref. m- 


3 Free Energy 



2r^ N (r)a( : r)f\r) _ 

—2r 2 N(j-)a'(r) + Ar 2 ^/N(r)a(r) 


—2rN(r)a(r)], (32) 


where U 2 denotes the area of the two dimensional trans¬ 
verse space. 

For the normal phase which is dual to the AdS-RN 
black brane solution 

H r ) = M (1 ^ —) > 

r 

3 9 9 

rv ‘J /" rr‘ /v» 

N(r) = r 2 (l-^)+b 2 fJ -l(l--), 

r A rh 

<?(r) = 1 , f(r) = 1 - ( 33 ) 

the temperature and free energy are 


In our model, there left one degree of freedom denoted 
by T for the solutions, once we have fixed the values of b 
and A. The model includes four different solutions, which 
are the AdS Reissner-Nordstrdm (RN) black brane solu¬ 
tion describing the normal conductor phase, the s-wave 
superconductivity solution, the p-wave superconductivity 
solution, as well as the s+p superconductor solution with 
both the expectation values for the scalar and vector or¬ 
der parameters non zero. These four kinds of solutions are 
all possible solutions with one free parameter T at fixed 
value of ( b,A), but the different solutions are in different 
branches. Thus at some fixed value of T, there may exist 
more than one solution, which means that the system may 
have more than one possible phase simultaneously. So we 
need to calculate the free energies of the different phases 
(solutions) to find out which is the most favored one. 

The free energy of the system is equal to the temper¬ 
ature T times the Euclidean on shell action of the bulk 
spacetime 

n = ts e . (30) 

In Se, we should consider both the bulk terms and some 
boundary terms. Remember that all the solutions are source 
free and we have fixed the chemical potential, which means 
we study the system in grand canonical ensemble. In this 
case, the action Se can be expressed as 


S E — 


- D^ a D^ a - m 2 & a ip a ) 


- 77 / d 3 x\/h[ K -). 


g J t — y oo 


(31) 


By substituting the solutions to the Euclidean action, the 
integrand can be written as a total derivative term, thus 
the on shell action and the free energy can be evaluated 
by the boundary value of the functions as 


T - —(3 - 2 k 9 q - r 3 1 b 2 u 2 r, 

~ 4tt (3 2r 2 ] ’ V 2 h 2 M h ’ 

while for the condensed phases, we have 

T = ^( 3 a h0 - ^ - b 2 m 2 a h0 & 2 h0 ) , 


(34) 


^n = -2 M b0 . 


(35) 


Using the above formulas, we can get the free ener¬ 
gies of all the solutions and compare these free energies 
at the same temperature and ( 6 , A) . Thus we can draw 
the condensate of the s-wave and p-wave orders for the 
most stable phase at any temperature with fixed b and 
A (Figures [lj [2| 3j |4j [5]). We can see what kind of phase 
transitions wouldoccur from these figures. At some values 
of b and A, the condensation behavior is very interesting, 
we show some typical cases in the next section. 


4 Condensation and phase transition 

We study the model with eight discrete values of the back 
reaction strength b from 0.1 to 0.8. In each value of 6 , 
we can fix the value of A to study the phase transitions 
of the system with a decreasing temperature. We should 
find all the possible solutions and calculate the relevant 
free energy to figure out the most stable phase at each 
temperature. Finally we can plot the condensation behav¬ 
ior of the s-wave and p-wave orders versus the temperature 
to show the phase transition at the fixed value of b and 

A. 

We can also change the value of A continuously to see 
what would happen to the stable phases and the phase 
transitions between them. With the information at differ¬ 
ent values of A, we can extract a phase diagram on the 
A — T plane at the fixed value of b. We will show the eight 
phase diagrams with different values of b in the next sec¬ 
tion. In this section we show some interesting and typical 
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Fig. 1. The condensation value of the s-wave and p-wave orders. The red line denotes the value for the p-wave order, and the 
blue line denotes the value for the s-wave order. Dotted lines are for unstable phases and solid lines for stable phases. 



T/n T/fi 


Fig. 2. The condensation value of the s-wave and p-wave orders. The red line denotes the value for the p-wave order, and the 
blue line denotes the value for the s-wave order. Dotted lines are for unstable phases and solid lines for stable phases. The right 
figure is an enlarged version for the rectangle region in the left figure. 


phase transition behaviors through the figure of conden¬ 
sation values versus temperature at fixed b and A. 

For each value of b, we vary the value of A to find 
all possible phase transitions. Below we pick some typical 
values of b and A 1 and draw figures for the condensation 
values of the s-wave and p-wave operators versus temper¬ 
ature in Figures [I] [2] [3] [3 [5] In these figures, we use the 
solid red line to show the condensation value of the p-wave 
operator, and the solid blue line to show the condensation 
value of the s-wave operator. We also use dotted lines to 
show the condensation values in unstable branches. The 
dotted red lines denote the condensation value of the p- 
wave order in the unstable section of the p-wave phases, 
and the dotted blue lines denote the value of the s-wave 
order in the unstable section of the s-wave phases. The 
s+p coexisting phase might also become unstable in our 
study. So in the unstable s+p phases, we use dotted or¬ 
ange line to denote the condensation value of the p-wave 
order and use the dotted green line to denote that of the 
s-wave order. In Figure[5]we also use vertical dashed black 
lines to denote critical temperatures of phase transitions 
between different condensed phases. 

In Figure^ the left plot is for the case b = 0.1 and A = 
1.9. In this figure, we can see a typical behavior of the s+p 
system. In this case the s+p coexisting phase connects the 
p-wave phase in higher temperatures and the s-wave phase 
in lower temperatures. This behavior is similar to that in 



Fig. 3. The condensation of the s-wave and p-wave orders. 
The red line denotes the value for the p-wave order, and the 
blue line denotes the value for the s-wave order. Dotted lines 
are for unstable phases and solid lines for stable phases. 

the probe limit [14j . but here the temperature region of 
the s+p coexisting phase is larger than that in the probe 
limit. 

The right plot of Figure [I] shows the case with b = 0.2 
and A = 1.93. We can see that the s+p coexisting phase 
emerges from a p-wave phase and is stable in the low tem¬ 
perature region (the condensate of the p-wave order does 
not vanish even in low temperature). In the s+p coexisting 
phase, the condensation value of p-wave order is reduced 
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when the value of the s-wave order increases from zero. 
This means that the s-wave order and the p-wave order 
still repel each other. 

Figure [2] shows the results with b = 0.3 and A = 1.885. 
The right plot is an enlarged version of the rectangle region 
in the left one. In this case, we find an interesting behav¬ 
ior for the s+p coexisting phase. The s+p coexisting phase 
begins from an s-wave phase at a higher temperature and 
finally goes back to the s-wave phase at a lower tempera¬ 
ture. We can see from the picture that the condensation 
value of the p-wave order forms a shape of the letter “n”, 
so we call this as an “n-type” condensation. In the n-type 
phase transition, the order parameter of the p-wave only 
gets non-zero values in a small region. Note that we can 
see from the right plot that the value of s-wave order in 
the s+p coexisting phase is reduced, compared to the case 
of the pure s-wave phase. 

In Figure [3j we show the results with b = 0.4 and 
A = 1.8544. In this figure, we can see that similar to the 
case in the right plot of Figure[l] the s+p coexisting phase 
is also stable in the low temperature region, but at this 
time, the s+p coexisting phase emerges from an s-wave 
phase instead. 

Figure [4] shows two similar cases. The left plot is for 
b = 0.5, A = 1.945 and the right one for b = 0.6, A = 
1.963. There are five regions of different phases in each of 
the plot, from the right to the left, they are the normal 
phase, the s-wave phase, the s+p coexisting phase, the 
p-wave phase and the second branch of s+p coexisting 
phase. The main difference between the two cases is that 
in the left plot, the phase transition from the s-wave phase 
to the s+p coexisting phase is second order, but in the 
right plot, this phase transition is first order. We can see 
from the two plots that the system goes into the p-wave 
phase from the s+p coexisting phase and goes back to the 
s+p coexisting phase at a lower temperature. The shape 
of the solid blue line forms a letter “u”, so we call this as 
a “u-type” condensation. These kinds of phase transitions 
are known as “reentrant” phase transitions in condensed 
matter physics, and have also been found in holographic 
study in Ref. ‘20» . 

Figure[5]is for the cases with b = 0.8 and three different 
values of A. In these plots, we use the vertical dashed 
black lines to denote the critical temperature of the phase 
transitions between different condensed phases. The first 
plot shows the case A = 2.1. We can see that in this 
case, there is no stable s+p phase and a first order phase 
transition occurs between the s-wave phase and the p-wave 
phase. However, an unstable s+p solution exists, and the 
condensation values of s-wave and p-wave orders in the 
unstable s+p phase are denoted by the dashed green and 
dashed orange lines. 

The second plot of Figure[5]presents the case with A = 
2.023. In this figure, we can see that there is also a hrst 
order phase transition between the s-wave phase and the 
p-wave phase. And similar to the case with A = 2.1, there 
is an unstable s+p solution connecting the s-wave phase 
and the p-wave phase. Besides this unstable s+p solution, 
there is another branch of stable s+p phase in the low 



Vi* 



Till 



Fig. 5. The condensation value of the s-wave and p-wave or¬ 
ders. The red line denotes the value for the p-wave order, and 
the blue line denotes the value for the s-wave order. Dotted 
lines are for unstable phases and solid lines for stable phases. 
The dotted orange line and dotted green line denote the con¬ 
densation value for the p-wave and s-wave orders in the unsta¬ 
ble part of the s+p solution. We use the vertical dashed black 
lines to mark the critical temperatures between the condensed 
phases clearly. 


temperature region. It is also interesting to notice that the 
two s+p phases overlap in some temperature region, which 
means one might find two different branches of solutions 
for the s+p coexisting phase at some temperatures. We 
should also notice that the condensation value of the p- 
wave order in the lower temperature branch of s+p phase 
differs from the value in the p-wave phase, although the 
two red lines (solid one for the s+p phase and dotted one 
for the p-wave phase) are very close to each other. 

The third plot in Figure [5] shows the case with A = 2. 
We can see that the p-wave phase in this case is always 
unstable while the s+p phase is partly unstable, and the 
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Fig. 4. The condensation value of the s-wave and p-wave orders. The red line denotes the value for the p-wave order, and 
the blue line denotes the value for the s-wave order. Dotted lines are for unstable phases and solid lines for stable phases. The 
dotted orange line and dotted green line in the right figure denote the condensation value for the p-wave and s-wave orders in 
the unstable s+p solution. 


phase transition from the s-wave phase to the s+p phase 
is first order. 

From the above figures we can see that this holographic 
system exhibits a rich phase structure. These different 
phase transition behaviors can be used to look for and 
test universal behaviors in holographic systems. In partic¬ 
ular, we have found that in this model, with some choice 
of the value of parameters, an s+p solution might be un¬ 
stable, which is quite different from the previous studies, 
where once an s+p solution appears, it is usually stable. 
In the next section we will construct the phase diagram 
of the system in terms of the dimension of the scalar or¬ 
der and the temperature. Before that, in the rest of this 
section, we will discuss the potential role of the unstable 
s+p solution in phase transitions between the s-wave and 
p-wave phases. 

The Swallow Tail for First Order Phase Transition 

Note that to get the above condensation behaviors, we 
have calculated the free energies for all the phases to con¬ 
firm the stability in each case. The critical temperatures 
for the first order phase transitions are also calculated 
from the free energy of the two related phases. From the 
free energy curves, we found that an s+p solution always 
exists when the two free energy curves for the s-wave phase 
and p-wave phase intersects. But this s+p solution is not 
always stable. We classify the three cases for the stability 
of this s+p solution as follows: 

1. Stable: In this case, there exist two second order phase 
transitions to and from the s+p phase, instead of a first 
order phase transition between the s-wave and p-wave 
phases. 

2. Totally Unstable: A first order phase transition be¬ 
tween the p-wave phase and the s-wave phase occurs. 

3. Partly Stable & Partly Unstalbe: A first order phase 
transition occurs between the s+p phase and the s- 
wave or p-wave phase. 

In our previous work Mi we have shown the free en¬ 
ergy curves for Case [T] where this s+p solution is stable. 


It is natural that in this case the s+p solution exists and 
is the most stable phase in its temperature region. But 
in the last two cases, especially in the second case, the 
s+p solution becomes totally unstable. Then we may ask 
what is the meaning for the existence of the unstable s+p 
solution? 

To answer the question, we draw the free energy curves 
in Figure [6] for the three condensation behaviors shown in 
Figure [5j The two plots in the first line are for the case 
b = 0.8, A = 2.1, the two plots in the second line are for 
the case b = 0.8, A = 2.023, and the last two figures are 
for the case b = 0.8, A = 2. In the plots on the left, we 
use the solid red line and solid blue line to denote the 
free energy of the p-wave phase and s-wave phase respec¬ 
tively. The green line denotes the free energy for the s+p 
phase(solution) and the dashed black line denotes that for 
the normal phase. We also draw the enlarged version for 
the last two cases to see which one is lower when the green 
and red lines are very close to each other in some region. 
We use the same data to draw the plots on the right side 
and those on the left side, and we use thick black curves 
to denote the swallow tail shape in first order phase tran¬ 
sitions in the plots on the right side. 

We can see that when the s+p solution near the free 
energy intersection point becomes partly or totally unsta¬ 
ble, a first order phase transition occurs. This first order 
phase transition can be from the s-wave phase to the p- 
wave phase as in the case when the s+p solution is totally 
unstable, or from the s-wave phase to the s+p phase as 
in the case when the s+p solution is partly unstable. In 
both the two cases, we can find a classical swallow tail in 
the free energy curves for the first order phase transition. 
To make this clear, we use thick black curves to mark the 
swallow tail shape in the plots on the right side of Fig¬ 
ure E Remember that we use the same data to draw the 
left and right figures on the same row. One can compare 
the left and right plots to see that although the unstable 
part of the s+p solution is not experienced in the phase 
transition, it is still quite necessary in forming the swal¬ 
low tail shape of the free energy curve. It seems that the 
existence of the unstable s+p solution is to bridge the s- 
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Fig. 6. The free energy of the different phases with b = 0.8. The first two figures are for the case A = 2.1, the third and forth 
figures are for the case A = 2.023 and the fifth and sixth ones for A = 2. In the three figures on the left side, the dashed 
black line denotes the free energy for the normal phase, the solid red line and solid blue line denotes the free energy of the 
p-wave phase and s-wave phase respectively, the solid green line denotes the free energy of the s+p phase (solution). We also use 
enlarged version in the third and fifth figure to show the detailed relation between the red and green lines when they are close 
to each other. In the three figures on the right side, we use thick black curves to show how the free energy curve of different 
phases form the shape of a swallow tail. 


wave phase to the p-wave phase and to form the standard 
swallow tail shape of free energy for the first order phase 
transitions. 

According to the previous results, the swallow tail shaped 
free energy curve seems always existing in any first order 
phase transition. If we assume the swallow tail shape of 
free energy curve as a universal property for the first or¬ 
der phase transition, we can explain the existence of s+p 
solution near the free energy intersection point as follows. 
When the free energy curves of the s-wave phase and the 
p-wave phase have an intersection point as in Figure [6j 
the s+p solution must exist. Otherwise a first order phase 
transition would occur between the s-wave phase and p- 
wave phase, and without the s+p solution, the free energy 
curve can not form a swallow tail shape, which contra¬ 


dicts with the universality we assumed. This contradiction 
could be solved by an s+p solution near the intersection 
point of free energy curves. This s+p solution could either 
be (partly) stable to avoid the first order phase transi¬ 
tion between the s-wave phase and p-wave phase, or be 
totally unstable but bridge the s-wave and p-wave phases 
and form the swallow tail shape of free energy curve for 
the first order phase transition between the s-wave and 
p-wave phases. We can also use the same logic to explain 
the existence of the s+d phase in the holographic s+d sys¬ 
tem studied in Ref. m- If our assumption for the swallow 
tail shape in first order phase transitions is true, it can be 
very useful in finding new solutions near some free energy 
intersection points in more general systems. 
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5 Phase diagrams 

In the previous section, we have shown the condensation 
behaviors for some interesting cases. If we get all the con¬ 
densation behaviors for different values of A at fixed b , 
we can summarize the information at different A to build 
a A — T phase diagram at this fixed value of b. In this 
section, we use the above method to give the A — T phase 
diagrams with eight different values of the back reaction 
parameter b from 0.1 to 0.8. With these phase diagrams, 
we can understand the rich phenomena shown in the previ¬ 
ous section much more clearly. And we can also compare 
the phase diagrams to study the effect of back reaction 
concretely. 

We show all the phase diagrams in Figure [7] and Fig¬ 
ure [8j respectively. In the eight phase diagrams with b 
from 0.1 to 0.8, we can see that all the diagrams have 
four different phases. These four phases are the normal 
phase, the s-wave phase, the p-wave phase and the s+p 
coexisting phase colored by light gray, light blue, light red 
and white, respectively. The thermodynamic stability of 
all these phases has been confirmed by our calculation of 
free energy. 

We also use different lines in the figures to show the 
sets of different phase transition critical points or emer¬ 
gent points of unstable phases. The solid lines divide the 
different phases, and the dashed lines near the zero tem¬ 
perature region are obtained by extrapolation from the 
solid lines with the same color. The dotted lines denote 
the points where an unstable solution emerges. We use 
different colors to mark the different critical lines. The sec¬ 
ond order s-wave and p-wave phase transitions from the 
normal phase have been well studied, we color the criti¬ 
cal lines black. When the back-reaction is strong enough, 
the p-wave phase transition becomes first order, so we use 
cyan line to denote the critical points. The s+p coexist¬ 
ing phase generally starts from the s-wave or the p-wave 
phase, so we use blue lines to denote the emergent points of 
the s+p coexisting phase on the s-wave phase, and use red 
curves to denote the emergent points of the s+p coexisting 
phase on the p-wave phase. When the red or the blue curve 
becomes dotted, the points denote the emergent points of 
unstable (part of) s+p solutions. In that case, a first order 
phase transition occurs. We use the solid yellow curves to 
denote the critical points of first order phase transitions 
between the s-wave phases and the p-wave phases, and use 
the solid green curves to denote the critical points of first 
order phase transitions between the s-wave phases and the 
s+p coexisting phases. 

In order to describe these phase diagrams more clearly, 
we define some special values of the dimension of the scalar 
order to mark vertical coordinates of some critical points. 
We list their definitions as follows. 

1. A p o'. The coordinate of the red line at zero tempera¬ 
ture. 

2. A s o: The coordinate of the blue line at zero tempera¬ 
ture. 

3. A t (A t+ , A t -): The coordinate of the intersection point 

of three or four solid lines. For the last three figures, 


there are two such intersection points, we denote its 
vertical coordinate value as A t+ and A t _, respectively, 
with A t+ > A t ,-. 

4. Aj: The coordinate of the contact point of the blue 
line and the green line. 

5. A pm i n , A sm i n : These two denote the coordinates for 
the lowest points on the red line and blue line respec¬ 
tively. But these two are only useful in some limited 
cases, for example, A pm i n in b = 0.5, 0.6 and A sm i n in 
6 = 0.3. 

Next we will discuss each phase diagram case by case. 
In each phase diagram, the vertical coordinate is the scal¬ 
ing dimension of the s-wave order parameter, and the hor¬ 
izontal coordinate is the temperature over the fixed value 
of chemical potential T/fi. In order to explain the phase 
diagrams more clearly, we can fix the operator dimension 
to a special value A = that means we choose a spe¬ 
cial horizontal line in the phase diagram, with the vertical 
coordinate value equal to A e , g .. We then force the system 
to be cooled slowly from a high temperature state in the 
normal phase region, so the point representing the phase 
of the system will go from the right side to the left along 
the line with A = A eg . We can see how many phase 
transitions would occur from the phase diagram, and we 
can also get the critical temperature as well as the order of 
each phase transition from the phase diagram. Notice that 
each condensation figure in the previous section shows the 
case of a horizontal line in one of the phase diagrams. 

Now let us begin from the first phase diagram with 
b = 0.1. The structure of this phase diagram is quite sim¬ 
ple, it is very similar to the phase diagram in the probe 
limit m- We can see that in this case Z\ p o > A s0 > A t . 
We can analyze the system in four regions of the param¬ 
eter A e , g ,. When A e , g , < A tl the system undergoes only 
one phase transition from the normal phase to the s-wave 
phase and stays in the s-wave phase in the low temper¬ 
ature region. When A t < A e , g , < Z\ s0 , there are three 
phase transitions. The system will first undergo a p-wave 
phase transition from the normal phase, then experience 
a second order phase transition from the p-wave phase to 
the s+p coexisting phase, and finally meet the final phase 
transition and change from the s+p coexisting phase to 
the s-wave phase. The left plot in Figure [l] shows the typ¬ 
ical condensation behavior of the system in this region. 
When Z\ s0 < A e , g . < A p0 , the system will also experi¬ 
ence the p-wave phase transition and the phase transi¬ 
tion from the p-wave phase to the s+p coexisting phase, 
but without the phase transition from the s+p coexisting 
phase to the s-wave phase, the system will finally stay in 
an s+p coexisting phase in the low temperature region. 
When A e . g . > A p0 , there is no s+p coexisting phase, and 
the the system will only undergo the p-wave phase transi¬ 
tion and keep in the p-wave phase in the low temperature 
region. 

The b = 0.2 phase diagram looks also very similar to 
the case in the probe limit. The only difference is that the 
region for the s+p coexisting phase is larger. There are 
also four regions for the parameter A, in these regions, 
the phase transition properties are qualitatively the same 
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Fig. 7. The A — T phase diagrams for b values from 0.1 to 0.4. In these figures, the areas colored by light gray, light blue, light 
red and white denote the normal phase, the s-wave phase, the p-wave phase and the s+p coexisting phase respectively. The 
black curves are made up by phase transition critical points between the s-wave or p-wave phases and the normal phase. The 
solid red curves describe the second order phase transition critical points between the p-wave phases and the s+p coexisting 
phases. The solid blue curves represent the second order phase transitions between the s+p coexisting phases and the s-wave 
phases. The dashed lines near the zero temperature regions are extrapolated from data in higher temperature region. 


as that in the case b = 0.1. We plot the condensation 
behavior of the system with A s o < Z\ e . g . < A p q in the 
right plot of Figure [lj 

In the b = 0.3 case, the phase diagram still looks simi¬ 
lar to the previous two cases. But it seems that the slope 
of the blue line changes to be positive in some region. We 
draw an enlarged version of the phase diagram to show 
the case clearly. From the enlarged version we can see 
that the blue line has a minimum value of A = A smin at 
a non zero temperature, below this temperature the blue 
line has a negative slope. We have four special values of 
A with A smin ^ ^sO ^ ^ ^pO) th.US WG llcLVG five 

regions for Z\ e . 5 .. The most interesting one is the region 
with Asr)tin ^e ■9- < Z\ s o, and the condensation behav¬ 
ior for this case is shown in Figure [2j In this region, there 
would be three phase transitions. The system first goes 
from the normal phase to the s-wave phase via an s-wave 
phase transition, then it goes from the s-wave phase to the 
s+p coexisting phase by the second phase transition. In 
the third phase transition, the system goes back to the s- 
wave phase from the s+p coexisting phase. In this case, the 
behavior of the p-wave order is rather interesting, the con¬ 
densation value of the p-wave order emerges from zero at 
some critical temperature. Rather than increasing rnono- 
tonically, it decreases after reaching its maximal value and 


finally vanishes when the system goes back to the s-wave 
phase. We plot the condensate of the s-wave and the p- 
wave orders for this kind of phase transition in Figure [2] 
and call this condensation as an “n-type” one. From the 
phase diagram we can see that the “n-type” condensation 
occurs because the blue line has a minimum value of A at 
a non-zero temperature. 

In the b = 0.4 case, the phase diagram is qualitatively 
the same as that in the case with b = 0.3, so the above 
analysis of the b = 0.3 phase diagram is still valid here. 
However, it should be pointed out that the lowest point 
of the blue line seems to be in the dashed part of the line. 
Because the dashed part is obtained by extrapolation, the 
lowest point needs to be confirmed by further numerical 
calculations. 

In the b = 0.5 phase diagram, the blue line always has 
a positive slope, but the red line has a minimum value 
A = A pm i n at a non-zero temperature in this case. Then 
there are four special values of A which will divide the 
phase diagram into five different regions. The interesting 
part is the region A pm i n < Z\ e . s . < A t . In this region, 
the s-wave phase transition occurs first, then the system 
goes into the s+p coexisting phase via a second phase 
transition. In the third and forth phase transitions, the 
system goes into the p-wave phase and then goes back 



















Zhang-Yu Nie et al.: Phase transitions in a holographic s+p model with back-reaction 


11 




T//i Till 


Fig. 8. The A — T phase diagrams for b values from 0.5 to 0.8. In these figures, the solid cyan line denotes the critical points of 
the first order phase transitions from the normal phase to the p-wave phase. The solid yellow lines represent the first order phase 
transition critical points from s-wave phases to p-wave phases. The solid green lines represents the first order phase transition 
critical points from s-wave phases to s+p coexisting phases. The dotted part of the red and blue lines are the emergent points 
of the unstable part of s+p phases. Notations of other lines are the same as introduced in Figure [7] 


to the s+p coexisting phase. We plot the condensation for 
this case in Figure |4j and we call this kind of condensation 
as “u-type”. We can see easily from the phase diagram 
that the “u-type” condensation can occur because the red 
line has a minimum value of A at a non-zero temperature. 

In the above five phase diagrams, there is always a 
quadruple intersection point connecting to all the four 
phases, and all the phase transitions are second order. In 
the next three phase diagrams, there does not exist any 
quadruple intersection point, but there are two triple in¬ 
tersection points. And we can see that in this case, some 
phase transitions become first order. 

In the phase diagram with b = 0.6, the p-wave phase 
transition from the normal phase is still second order, but 
in the region Aj < Z\ e . s . < A t -, the phase transition 
from the s-wave phase to the s+p coexisting phase be¬ 
comes first order. When A t - < A e , g , < A t+ . the s+p 
coexisting phase connecting the s-wave phase and p-wave 
phase still exists, but it always has a higher free energy 
than the s-wave and p-wave phases. Thus a first order 
phase transition occurs between the s-wave phase and the 
p-wave phase. But as we have discussed in the previous 
section, the unstable s+p coexisting phase is still essential 
to form the swallow tail shape of the free energy curve for 
the first order phase transition. In this phase diagram, the 
red line still has a minimum point as in the b = 0.5 phase 



Fig. 9. The enlarged version of the A — T phase diagrams for 
the case with b = 0.6. 


diagram, then the system would also go into the p-wave 
phase from an s+p coexisting phase and goes back to the 
s+p coexisting phase later when A pm i n < A e g . < A t -. 
We plot the condensation behavior for this case in the 
right plot of Figure [4] We draw an enlarged version of this 
phase diagram in Figure [9] to show the details near the 
triple intersection points. 

When b = 0.7, the p-wave phase transition from the 
normal phase becomes first order, so the vertical line sepa- 
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rating the p-wave phase and normal phase is marked cyan. 
In this case, the s+p phase transition from the s-wave 
phase is still first order in the region Aj < A eg , < At-, 
and the s+p coexisting phase connecting the s-wave and 
p-wave phases is always unstable in the region A t - < 
Ae.g . < A{+. 

The b = 0.8 phase diagram is very similar to the one 
with b = 0.7, but notice that the green line seems not 
to intersect with the blue line, which means the phase 
transition from the s+p coexisting phase to the s-wave 
phase will be alway first order. We show the condensate 
behaviors in Figure [5] for three typical regions for A eg . 
In Figure [5] the first plot is for A p0 < A e , g _ < A t+ , the 
second plot is for the case A t - < A e _ g , < A p0 , and the 
third one is for Zi s o < A e . g . < At-. We also plot the free 
energy for the three cases in Figure [6] 

From the above analysis, we can see that the qualita¬ 
tive condensation behaviors at fix b and A can be read off 
easily along a horizontal line in each phase diagram. The 
different shape of the curves separating the four phases 
leads to different phase transition behaviors. For example, 
the “n-type” behavior in Figure [2] is due to the existence 
of the minimum of the blue line in the b = 0.3 phase dia¬ 
gram, while the “u-type” behavior in Figure[4]is due to the 
existence of the minimum of the red lines in the b = 0.5 
and b = 0.6 phase diagrams. 

In order to better understand the above phase transi¬ 
tion behaviors, we try to find out the possible mechanism 
behind them. Because the phases with individual s-wave 
or p-wave order have been well studied, the key issue is 
to find the underlying law for the emergence of s+p co¬ 
existing phase. However, our study is based on numeral 
solutions, thus we can only give some heuristic argument 
as follows. 

In the holographic superconductor model with single 
s-wave or p-wave order parameter, the condensate of the 
order begins when the electric field near the black brane 
horizon is large enough to support formation of the charged 
hair. If we lower the temperature with fixed value of the 
chemical potential in the normal phase, the electric field 
near the horizon is increasing. Therefore, after the strength 
of the electric filed is larger than some critical value, the 
charged order will condense. In the case of our s+p model, 
the s-wave and p-wave order generally have different crit¬ 
ical temperature, therefore the one with higher critical 
temperature will condense at first. When the first order 
begins to condense, the condensed order will deplete some 
charge from the horizon. As a result, the strength of the 
electric field near horizon is reduced to a smaller value 
than that of the normal phase at the same temperature. 
Therefore the condensate of the second order in presence 
of the first order should occur at a lower temperature than 
the condensate of the second order on the normal phase 
(without the first order). We can also confirm this from 
the phase diagrams in this section. 

Another important quantity in the study of the co¬ 
existing phase is the free energy. We find that the s+p 
phase always exists once the free energy curves for the 
s-wave phase and p-wave phase intersect. It might be pos¬ 


sible to get a qualitative formula to predict the existence of 
the s+p phase, by using some information in free energy. 
However, this need to be further investigated in future. 

We can also see the influence of the back reaction from 
the phase diagrams. The most obvious one is that the 
region for the s+p coexisting phase is enlarged by the 
back reaction when the back reaction strength b is not too 
large. This is in accordance with the intuition that the 
back reaction on the metric adds an additional attraction 
between the s-wave and p-wave orders. This also happens 
in the s+s system in Ref. m- 

However, we here have found that for a very large value 
of b, the region for the s+p coexisting phase is reduced. 
This is mainly caused by the low critical temperatures 
of the s-wave and p-wave phase transitions in the strong 
back-reaction cases. We should also notice that the s+p 
solution near the triple intersection points becomes unsta¬ 
ble in the strong back reaction case, this also makes the 
region for the s+p coexisting phase smaller. 

In the parameter region we have tested, the s-wave 
phase transition from the normal phase is always second 
order, but the p-wave phase transition becomes first order 
when the back reaction is large enough. The critical value 
of the back reaction at which the p-wave phase transition 
becomes first order can be read off from Ref. [.'12] . and 
the critical value is b c = 0.62. In this work, we see that 
the phase transition from the s-wave phase to the s+p 
coexisting phase can also be first order when the back 
reaction is large. And even when both the p-wave and 
s-wave phase transitions are still second order, the s+p 
coexisting phase transition can be first order, as exhibited 
in the phase diagram of b = 0.6. 

6 An alternative setup of the holographic 
s+P model 

Before we conclude our work, we give another setup for 
a holographic s+p model in this section. From this new 
setup, we can get the same equations of motion as those 
presented in Sec. [2] Thus we can get the same conden¬ 
sation behavior and the same phase diagrams as in the 
previous sections. 

It is well known that the origin of superconductivity 
involves the formation of a quantum condensate state by 
pairing conduction electrons. Nevertheless, from a symme¬ 
try point of view, the phenomenon of superconductivity 
can be thought of as the spontaneously breaking of U(l) 
gauge symmetry. In the setup discussed above, the p-wave 
and s-wave order parameters are dual to the gauge field 
pointing in the third direction inside the SU(2), i.e., A 3 
and the third component of the scalar triplet 1 F 3 , respec¬ 
tively. By fixing a gauge, we choose the SU(2) gauge field 
in the first direction as the U(l) gauge symmetry which 
we call U(l)i. A non-vanishing value of the time com¬ 
ponent of 1/(1) 1 gauge field induces a chemical potential 
on the dual theory. The non-trivial profile of the gauge 
field A 3 would induce a vacuum expectation value of the 
dual operator with no source and thus breaks the C7(l)i 
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symmetry spontaneously. Similar story is also true for the 
scalar field F 3 . 

However, there are still some limitations in the setup 
discussed above. One can change the scaling dimension of 
the dual scalar operator for the s-wave order by adjusting 
the mass square of F a . But dual to the constraint from the 
SU(2) gauge symmetry, the dimension of the p-wave order 
has to be fixed to be A p = d — 1 where d, is the spatial 
dimension in the bulk. Furthermore, the SU(2) gauge in¬ 
variance also makes the “charge” of the p-wave order equal 
to the one of the s-wave order. Therefore in our previous 
setup, we only have two free parameters b and A s = A 
at hand. In principle, the dimension and the charge for 
the p-wave order can also have different values, so we also 
need to find the way to vary these parameters in the holo¬ 
graphic setup. One way is to use the holographic p-wave 
model [30] in terms of a charged complex vector field. In 
this new setup, we can use a complex scalar and a complex 
vector minimally coupled with a 17(1) gauge field as the 
matter part of the holographic action. The matter part of 
the action can be written as 

S’m=J d d+ 1 xV^(-\F^F^ - 

-D^FD^F - m 1 2 s FF t). (36) 

where the superscript “t” means complex conjugate, F pv = 
S7 p A l , — „A p is the field strength for the U(l) gauge field. 
Here F is now a complex scalar charged under the U(l) 
gauge field with D p = S7 p — iq s A where q s is the U(l) 
charge of F. p p is a complex vector field charged under 
the same U(l) field but with a different charge q p . The 
field strength of p p is p pl , = D p p v — D u p p with covari¬ 
ant derivative D p = — iq p A p . The same equations of 

motion from the previous model under the ansatz ([9]) can 
also be obtained in the particular case with q p = q s = g c 
and m p = 0 in this new model under the following ansatz 

F = F(r), A t = 4>(r), p x = F x (r), (37) 

with all other matter field components being set to zero. 
Thus all results in the previous sections can also be rec¬ 
ognized as the q s = q p = g c and m p = 0 case in the new 
s+p model in this section. 

It is worth pointing out that although those two s+p 
models show same aspects in particular cases, they are 
distinct from each other. For example, the transport coef¬ 
ficients might exhibit essential differences [52] . And with 
this new holographic s+p model, we can consider the effect 
of different charges of s-wave and p-wave orders. Further¬ 
more, we can also introduce the mass term of p p to deform 
the scaling dimension of the p-wave vector operator. So we 
can study the s+p system holographically in a more com¬ 
plete parameter space. It has been shown that depending 
on q p and m p , only the p-wave part of action ( |36| ) can 
exhibit a rich phase structure mm- Thus the new 
model ( |36[ ) is expected to show more fruitful phase be¬ 
haviors |43L144] in a larger parameter space. This will be 
left for further study. 


7 Conclusions and discussions 

In this paper, we have continued our study on the holo¬ 
graphic s+p model with a scalar triplet charged under 
an SU(2) gauge field in the bulk [14] . by considering the 
back reaction of the matter fields on the background black 
brane geometry. We have shown that the model exhibits 
a rich phase structure depending on the model param¬ 
eters. Taking eight different values of the back reaction 
parameter b from 0.1 to 0.8, we have plotted the conden¬ 
sate behaviors of the s-wave and p-wave orders and built 
corresponding phase diagrams in terms of the dimension 
of the s-wave order and the temperature of the system. 
We have also given another setup for a holographic s+p 
model, which can lead to the same equation of motion as 
the model proposed in [14] with special value of model 
parameters. 

In this holographic s+p superconductor model, we have 
discovered some new condensate behaviors such as the 
“n-type” and “u-type”, where the condensations for the 
p-wave or s-wave orders are non-monotonic. These inter¬ 
esting behaviors can only occur when the blue or red curve 
in the phase diagram has a minimum at non-zero temper¬ 
ature. These non-monotonic condensation behaviors are 
also known as “reentrant’Qphase transitions in condensed 
matter physics, and would be very useful in future stud¬ 
ies. For example, we can solve the time evolution problem 
of a quenched initial state in the “n-type” to study the 
non-equilibrium properties in systems with multi-orders 
and compare the results to that in the system with only 
the s-wave order. 

We also have found that the s+p solution might be 
totally unstable when the back reaction is strong enough, 
and a first order phase transition between the s-wave and 
p-wave phases occurs at that time. Although the s+p so¬ 
lution in that case is totally unstable, it is still essential 
in forming the swallow tail in free energy curves for the 
first order phase transition. According to this property, if 
we assume the swallow tail of free energy curve as a uni¬ 
versal feature for the first order phase transitions, we can 
explain why the s+p solution can always exist when the 
free energy curves for the s-wave and p-wave phases have 
an intersection point. 

The eight phase diagrams give a unified description for 
the various condensation behaviors in Sec. [2| Each phase 
condensation behavior at fixed A and b can be read from 
a horizontal line in the phase diagram. We can also get a 
concrete understanding on the effect of the back reaction 
on our holographic s+p model from the resulting phase di¬ 
agrams. We can see that the back reaction of the matter 
sector on the black brane background generally enlarges 
the region for the s+p phases when the back reaction is 
not too large, which is similar to the case of the s+s sys¬ 
tem m- But in the very strong back reaction case, the 
region for the s+p phase becomes smaller. This is mainly 
caused by the decreasing of the critical temperatures for 
the s-wave and p-wave phase transitions when the back 

1 We would like to thank Professor Jan Zaanen for reminding 

us this issue. 
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reaction becomes strong. The fact that the s+p phase be¬ 
comes unstable near the triple intersection points also con¬ 
tributes to the decreasing of the region for the s+p phase 
in the phase diagram. 

In the strong back reaction case, the p-wave phase 
transition from the normal phase becomes first order. With 
the back reaction getting stronger, there is a growing pa¬ 
rameter space where the transition from the s-wave phase 
to the s+p phase is first order. The first order phase tran¬ 
sition from the s-wave phase to the s+p phase even exists 
when both the s-wave and p-wave phase transitions are 
still second order. 

There are many further directions to explore based on 
this work. With the new setup proposed in section [6j we 
can study the s+p system more completely with various 
charges and dimensions of the p-wave order. Since our 
s+p system has shown many new phase transition behav¬ 
iors, one can use them to test some universal properties 
in holographic systems. In addition, these new s+p mod¬ 
els provide the setups in which we can study time de¬ 
pendent or inhomogeneous solutions in order to see the 
non-equilibrium physics or the lattice effects. Finally if 
the swallow tail shape of the free energy curve is indeed 
a universal property for the first order phase transitions, 
it can be used to predict new phases(solutions) near the 
intersection point of two different phases. 
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